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' Abstract 

O 

X ; A search for lepton-flavor-violating interactions ep — > fiX and ep — > rX" has 

V | been performed with the ZEUS detector using the entire HERA I data sample, 

^| corresponding to an integrated luminosity of 130 ptT 1 . The data were taken 

at center-of-mass energies, \/s, of 300 and 318 GeV. No evidence of lepton- 
flavor violation was found, and constraints were derived on leptoquarks (LQs) 
that could mediate such interactions. For LQ masses below ^/s, limits were set 
on A egi y / ^, where \ eqi is the coupling of the LQ to an electron and a first- 
generation quark g 1; and (3e q is the branching ratio of the LQ to the final-state 
lepton i (yU or r) and a quark q. For LQ masses much larger than y^i, limits were 
set on the four-fermion interaction term X e q a ^eq fJ / '^lq f° r LQs that couple to an 
electron and a quark q a and to a lepton £ and a quark qp, where a and (3 are 
quark generation indices. Some of the limits are also applicable to lepton-flavor- 
violating processes mediated by squarks in i?-Parity-violating supersymmetric 
models. In some cases, especially when a higher-generation quark is involved 
and for the process ep — > tX, the ZEUS limits are the most stringent to date. 
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1 Introduction 



The recent observations of neutrino oscillations [1,2] have shown that lepton- flavor vio- 
lation (LFV) does occur in the neutrino sector. The LFV induced in the charged-lepton 
sector due to neutrino oscillations cannot be measured at existing colliders due to the low 
expected rate [3]. However, there are many extensions of the Standard Model (SM) such 
as grand unified theories (GUT) [4], supersymmetry (SUSY) [5], compositeness [6] and 
technicolor [7] that predict possible e — > fi or e — > r transitions at detectable rates. 

In many theories, LFV occurs only in the presence of a particular quark generation. At 
the HERA ep collider, lepton-flavor-violating interactions can be observed in the reaction 
ep — > £X, where I is a \x or r. The presence of such processes, which can be detected 
almost without background, would clearly be a signal of physics beyond the Standard 
Model. This search is sensitive to all quark generations for LFV occurring between e 
and fj, or r. Strong constraints on LFV also arise from measurements of rare lepton and 
meson decay, muon-electron conversion on nuclei, etc. [8]; nevertheless, HERA generally 
has a competitive sensitivity, and better sensitivity in the case of e - r transition when a 
second- or third-generation quark is involved. 

In this search, no evidence for LFV was found. The Buchmiiller-Ruckl-Wyler (BRW) 
leptoquark (LQ) model [9] and supersymmetry with i?-Parity violation are used to set 
limits from the search. Leptoquarks are bosons that carry both leptonic (L) and baryonic 
(B) numbers and have lepton-quark Yukawa couplings. Their fermionic number (F = 
3B + L) can be F = or \F\ = 2. Such bosons arise naturally in unified theories that 
arrange quarks and leptons in common multiplets. A LQ that couples both to electrons 
and to higher-generation leptons would induce LFV in ep collisions through the s- and 
w-channel processes shown in Fig. 1. The same processes can also be mediated by squarks, 
the supersymmetric partners of quarks, in SUSY theories that violate i?-Parity [10]. A 
detailed description of the considered phenomenological scenarios and of the cross section 
assumptions used in this paper is given in a previous publication [11]. 

Searches for LFV have been previously made at HERA [11,12]. This analysis is based on 
the entire HERA I sample collected by ZEUS in the years 1994-2000, corresponding to 
an integrated luminosity of 130pb _1 . These results supersede previous results published 
by ZEUS [11, 13], based on a sub-sample of the present data. 

2 The experimental conditions 

A detailed description of the ZEUS detector can be found elsewhere [14]. In this section 
a brief outline of the main components used in this analysis is given: the central tracking 
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detector (CTD) [15], the uranium-scintillator calorimeter (CAL) [16] and the forward 
muon detector (FMUON) [14]. 

The CTD, which is immersed in a magnetic field of 1.43 T provided by a superconducting 
solenoid, consists of 72 cylindrical drift chamber layers, organized in 9 superlayers covering 
the polar-angle 1 region 15° < 9 < 164°. The transverse-momentum resolution for full- 
length tracks is a{p T )/p T = 0.0058p T © 0.0065 © 0.0014/p T , with p T in GeV. The CTD 
was used to reconstruct tracks of isolated muons and charged r-decay products. It was 
also used to reconstruct the interaction vertex with a typical resolution of 4 mm (1mm) 
in the Z (X and Y) coordinate. 

The high-resolution uranium-scintillator calorimeter consists of three parts: the forward 
(FCAL), the barrel (BCAL) and the rear (RCAL) calorimeters. Each part is subdi- 
vided transversely into towers and longitudinally into one electromagnetic section (EMC) 
and either one (in RCAL) or two (in BCAL and FCAL) hadronic sections (HAC). The 
smallest subdivision of the calorimeter is called a cell. The CAL energy resolutions, 
as measured under test-beam conditions, are a(E)/E = 0.18/\/i? for electrons and 
a(E)/E = 0.35/ VE for hadrons (E in GeV). 

The FMUON detector, located between Z = 5 m and Z — 10 m, consists of 6 planes of 
streamer tubes and 4 planes of drift chambers. The magnetic field of 1.6 T produced by 
two iron toroids placed at about 9 m from the interaction point and the magnetic field of 
the iron yoke (1.4 T) placed around the CAL enable the muon-momentum measurements 
to be made. The use of FMUON extends the acceptance for high-momentum muon tracks 
in the polar-angle region 8° < 6 < 20°. 

The luminosity was measured using the process ep — > &yp. The small-angle photons were 
measured by the luminosity detector [17], a lead-scintillator calorimeter placed in the 
HERA tunnel at Z = -107 m. 

2.1 Kinematic quantities 

The total four-momentum in the CAL (E, Px, Py, Pz) is defined as: 

Cy j Ei, 22 E i sin #i cos 4>h ^2 Ei sin ^ sin z2 Ei cos 

i i i i 

where Ei is the energy measured in the i th calorimeter cell. The angular coordinates 9i 
and 0j of the i th cell are measured with respect to the reconstructed event vertex. The 

1 The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction" , and the X axis pointing left towards 
the center of HERA. The coordinate origin is at the nominal interaction point. 
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absolute value of the missing transverse momentum, P t , is given by a/ P\ + P Y > while the 
transverse energy, E t , is defined as ^^sin^j. 

Another relevant quantity used in this analysis is E — Pz = ^2iEi(l — cos9i). In the 
initial state, E — Pz = 2E e , where E e is the electron beam energy of 27.5 GeV. If only 
particles in the very forward direction (proton beam), which give negligible contribution 
to this variable, are lost, as in NC DIS events, E — Pz ~ 55 GeV is measured in the final 
state. 

Jets are reconstructed using the kx cluster algorithm [18] in the inclusive mode [19]; only 
jets with transverse momentum greater than 4 GeV are considered. 

3 Data samples and Monte Carlo simulation 

The data used in this analysis were collected in the years 1994-2000. The total integrated 
luminosity was 112.8 ± 2.2pb _1 with e + p collisions at the center-of-mass energy of 300 
and 318 GeV and 16.7 ± 0.3 pb" 1 with e~p collisions at 318 GeV. 

In the absence of a signal, limits were placed on LFV coupling strengths. The search is 
sensitive to any process with a final-state topology where the scattered electron of the ep 
neutral current (NC) deep inelastic scattering (DIS) is replaced with a fi or a r. However, 
for the purpose of limit setting, the signal was taken to be the LFV processes mediated 
by scalar or vector LQs of any mass. These were simulated by the Monte Carlo (MC) 
generator Lqgenep 1.0 [20], which is based on the BRW model. The simulation of the 
hadronization and particle decays was performed using Pythia 6.1 [21]. 

Various MC samples were used to study the Standard Model background. Charged cur- 
rent (CC) and NC DIS events were simulated using Djangoh 1.1 [22], an interface to 
the program Heracles 4.6.1 [23] and Lepto 6.5.1 [24]; Herwig 6.1 [25] was used for 
photoproduction background simulation while lepton pair production was simulated with 
Grape 1.1 [26]. 

4 e — \i transition 

The characteristic of such events is an isolated muon with high transverse momentum, 
which is balanced by that of a jet in the transverse plane. An apparent missing transverse 
momentum, measured by the calorimeter, due to the penetrating muon is used for event 
selection. Further requirements were applied, as described below, to identify charged 
particles as muons. 
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4.1 Muon identification 



The muon identification comprises two different methods, in two different angular regions, 
for the final-state \i candidate. The first was used in the polar-angle range 15° < 9 < 164° 
and required that the following conditions were satisfied: 

• a CTD track pointing to the vertex with transverse momentum above 5 GeV matching 
a calorimeter deposit compatible with a minimum-ionizing particle; 

• D trk > 0.5 and D jet > 1 where D t rk (-Djct) is the distance in the r\ — <fi plane between the 
track associated with the candidate muon and the closest track (jet) to the candidate; 

• candidate muons in the polar-angle region 115° < 9 < 130° were excluded to eliminate 
background from electrons that lose much of their energy in the dead material at the 
transition between BCAL and RCAL. 

The second method was used for very forward muons (8° < 9 < 20°) and required a 
reconstructed track in the FMUON detector with hits in at least 5 detector planes. 

4.2 Preselection 

The trigger used in this analysis was identical to that used in CC DIS measurement 
described in detail elsewhere [27]. It was based on a cut on f t with a threshold lower 
than that used in the offline analysis. After applying timing and other cuts to reject 
background due to non-ep collisions (cosmics and beam-gas interactions), the following 
preselection requirements were imposed: 

• a reconstructed vertex with |Z v t x | < 50 cm; 

• f t > 15 GeV; 

• no electron 2 candidate with energy larger than 10 GeV [28]; this cut was used to 
suppress NC DIS processes in a region of potentially high background and negligible 
anticipated signal; 

• an isolated-muon candidate in the direction of the f t (A0 < 20°, where A<fr is the 
difference between the azimuthal angles of the candidate muon and of the ft vector). 

After the preselection, the sample contained 20 data events, while 25.9±1.1 were expected 
from SM MC, mainly from QED di-muon processes (ep — > fi + fi~X). 

2 Throughout this paper, "electron" is used generically to refer to e + as well as e~. 
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4.3 Final selection 

The cuts for the final selection were designed optimizing the sensitivity using signal and 
background simulations [29]. For this purpose a scalar LQ with a mass of 600 GeV, 
coupling to second-generation quarks, was taken as signal. Monte Carlo studies showed 
that this procedure results in good sensitivity for the whole range of LQ masses considered 
here. The following cuts were applied: 

• f > 20 GeV; 

• ft/VEt > 3V GeV; this cut was chosen to reject high-i? 4 background events, where 
the small apparent ft can arise from the finite energy-measurement resolution; 

• E — P z + > 45 GeV, where A M = f(l — cos # M ) / sin 6^, 6^ being the polar angle of 
the candidate muon; the quantity A M represents the contribution to E — Pz carried 
by the muon, assuming that the transverse momentum of the muon is ft- 

Figure 2 shows the comparisons between data and MC expectations before the final selec- 
tion. No event satisfied the final cuts, while 0.87±0.15 were expected from the simulation 
of the SM background. 

For LFV events mediated by resonant production of a leptoquark, the selection efficiency 
varied with the LQ mass, ranging from 39% to 54% for scalar LQs and from 47% to 62% 
for vector LQs with mass between 140 and 300 GeV. 

For leptoquarks with mass much greater than the center-of-mass energy the efficiency is 
typically lower than that for resonant LQs, because of the softer Bjorken-x distributions 
of the initial-state quarks. In this case the efficiencies were almost independent of the LQ 
mass but depended on the generation of the initial-state quark. Sea quarks, with softer 
Bjorken-x distribution than valence quarks, result in a lower momentum of the final-state 
lepton, leading to a lower signal efficiency Overall, the selection efficiency for high-mass 
LQs was in the range 20-45%. 

5 e — r transition 

Lepton-flavor-violating events leading to a final-state r are characterized by a high- 
momentum isolated r balanced by a jet in the transverse plane. Since the r decays 
close to the interaction vertex, only its decay products are visible in the detector. Due 
to the presence of at least one neutrino in all r-decay channels, a high value of ft is 
expected. Therefore, for all the channels, the CC DIS trigger (as described in Section 4.2 
for the muon channel) was used together with the following common preselection: 

• f > 15 GeV; 
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• a reconstructed vertex with \Z vtx \ < 50 cm. 

5.1 Leptonic r decays 

For r leptons decaying into unions (r — > [w^v^), the same selection cuts as described in 
Section 4.3 were applied, since the event topology is very similar to that of LFV with 
e — > \i transitions. 

For the r — > ev e v T channel, the final state is characterized by a high-energy isolated 
electron in the ft direction; the following cuts were applied after the preselection: 

• 20 < E - P z < 52 GeV; 

• total energy deposit in RCAL less than 7 GeV; 

• p t /VE t > 2.5/GeV; 

• an electron with energy larger than 20 GeV in the polar- angle region 8° < 6 < 125° 
and in the f t direction (A0 < 20°); 

• a jet with a transverse momentum above 25 GeV, back-to-back with respect to the 
electron (A(p e ~^ ct > 160°) where A<f) e ~i Qt is the difference between the azimuthal angles 
of the jet and of the electron. 

No event was found in data, while 0.43 ± 0.08 were expected from SM MC. 

5.2 Hadronic r decays 

The r lepton, because of its small mass, typically decays with only one or three charged 
tracks with limited transverse spread. Since jets coming from hadronic r decays must be 
separated from a large background of QCD jets, a r finder was employed to distinguish 
the r jets from the quark- and gluon-induced jets. The algorithm exploits the fact that 
high-energy QCD jets usually have higher multiplicity and a larger internal transverse 
momentum than those for the decay products of the r. 

5.2.1 Tau identification 

A technique for r identification [30] was developed for a previous study [31] in which 
a small number of isolated-r events were found in the data set identical to that used 
here. The longitudinally invariant kx cluster algorithm was used to identify jets. The 
jet shape was characterized by the following six observables [31]: the first (-R mea n) and 
the second (-R rm s) moment of the radial extension of the jet-energy deposition; the first 



(3 



(-^mean) and the second (L rms ) moment of the energy deposition in the direction along the 
jet axis; the number of subjets (iV su bj) within the jet resolved with a resolution criterion 
y cut of 5 ■ 1CT 4 [32,33]; the mass (Mj et ) of the jet calculated from the calorimeter cells 
associated with the jet. In order to separate the signal from the background, the six 
variables were combined into a discriminant "D, given, for any point in the phase space 
x(- log(E mcan ), - log(i? rms ), - log(l - L mcan ), - log(L rms ), N suh j, M jct ), by: 

V(x) = Psig ^ , 

where p S i g and pbkg are the density functions of the signal and the background, respectively 
Such densities, sampled using MC simulations, were calculated using a method based on 
range searching [34]. Lepton-flavor- violating events in which the final-state r decays into 
hadrons and a neutrino were used to simulate the signal. The background simulation was 
based on CC DIS MC events. For any given jet with phase space coordinates x, the signal 
and the background densities were evaluated from the number of corresponding simulated 
signal and background jets in a 6-dimensional box of fixed size centered around x. The 
r signal tends to have a large discriminant value (T> 1) while the CC DIS background 
has a low discriminant value (T> — > 0). 



5.2.2 Preselection 

The following cuts were applied for the preselection of the hadronic r decay channel: 

• no electron candidate with energy larger than 10 GeV; 

• Et > 45 GeV; 

• 15 < E - P z < 60 GeV; 

• total energy deposit in RCAL less than 7 GeV; 

• a r-jet candidate as described below. 

The r-jet candidate was required to have a transverse momentum greater than 15 GeV, to 
be within the CTD acceptance (15° < 9 < 164°) and to have between one and three tracks 
pointing to the CAL energy deposit associated with the jet. Events with jets in the region 
between FCAL and BCAL (36° < 9 < 42°) were removed. In order to reject electrons 
from NC events, a cut of 0.95 was applied to the electromagnetic energy fraction of the 
jet (/emc)- In addition the jet was required to satisfy the condition /lt + /emc < 1-6, 
where /lt (the leading-track fraction) was defined as the ratio between the momentum of 
the most energetic track in the jet and the jet energy. The quantity /lt + /emc is close 
to 2 for electrons, the main source of background that this cut is designed to reject. 
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Figure 3 shows the comparison, after the preselection, between data and MC for the 
jet discriminant variables. Figure 4 compares the discriminant and the A0 distributions. 
Here, A0 is the azimuthal angle between the candidate r-jet axis and the f t vector. After 
the hadronic preselection, 119 events were found in the data, while 131 ±4 were expected 
from SM processes, mainly from CC DIS. The data distributions in Fig. 4 generally 
conform to those expected from SM backgrounds. 

5.2.3 Final selection 

For the final selection, the following additional cuts were applied to the events in Fig. 4: 

• V > 0.9; 

• the r-jet candidate was required to be aligned in azimuth with the direction of the f t 
(A0 < 20°). 

The discriminant cut was tuned to optimize the separation power, S = e sig • \fR (where 
e sig is the signal efficiency and R = 1/fbg is the background rejection), for a scalar LQ 
with a mass of 240 GeV [29]. In Fig. 5, the A0 distribution of the 8 events with T> > 0.9 
is shown compared to the SM expectation (10.2 ± 0.9 events). 

After imposing the final cut on A<f), no data events remained in the hadronic decay 
channel, while 1.1 ± 0.5 were expected from MC. 

5.3 Summary on e — > r search 

No candidate was found in the data for any of the three r-decay channels, while 2.3 ± 0.5 
were predicted by Standard Model simulations. 

The combined selection efficiency for low-mass (Mlq < \/s) scalar (vector) LQs was in 
the range of 22 - 29% (23 - 34%), while for high-mass (M LQ > y/s) LQs it was 4-20%. 
As is the case for the e — > \x transition discussed above, the significant efficiency drop for 
high-mass LQs is due to the softer Bjorken-x distribution of the initial state quarks. 

6 Results 

Since no evidence of lepton-flavor-violating interactions was found, limits at 95% C.L. 
were set - using a Bayesian approach [35] that assumes a flat prior for the signal cross 
section - on the processes ep — > fiX and ep — > tX mediated by a leptoquark. 
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In the low-mass case, limits on the cross section were converted, using the narrow-width 
approximation, into limits on A egi y / ^, where X egi is the coupling between the lepto- 
quark, the electron and a first-generation quark, while j3^ q is the branching ratio of the 
leptoquark into a lepton £ and a quark (u, d, s, c, b). For high- mass leptoquarks, the cross- 
section limits were converted, using the contact-interaction approximation, into limits on 
Aeg^A^/MLQ, where a and (3 are quark generation indices. The cross sections were eval- 
uated using the CTEQ5 [36] parton densities, taking into account the QED initial-state 
radiation, and, for low-mass scalar leptoquarks, NLO QCD corrections. 

6.1 Systematic uncertainties 

The following sources of systematic uncertainties are dominant: 

• the calorimeter energy-scale uncertainty (2%). The resulting variation in the signal 
efficiency for the muon (r) channel is less than 1% (3%) for low-mass leptoquarks and 
less than 5% for high-mass leptoquarks; 

• the luminosity uncertainty: 1.5% for the 1994-97 e + p data, 1.8% for the 1998-99 e~p 
data and 2.2% for the 1999-2000 e + p data; 

• Systematics related to the parton-density functions (PDF) have been calculated us- 
ing the 40 eigenvector sets, provided by CTEQ 6.1 [37], that characterize the PDF 
uncertanties. This contributes to the dominant uncertainty for low-mass leptoquarks, 
especially when a d quark is involved and the LQ mass approaches the HERA kine- 
matic limit. The effect of this uncertainty on the LQ limits is given in more detail 
elsewhere [29]. 

The uncertainties related to muon and tau identification were evaluated following the 
methods described elsewhere [31,38] and were found to be small. The systematic uncer- 
tainties have been included in the limit calculation assuming a Gaussian distribution for 
their probability densities. For low-mass LQs, the effect of the inclusion of systematic 
uncertainties is the largest at the highest masses and the limit on the coupling increases 
by less than 7% at 250 GeV. The effect is very small for high- mass LQs (below 1%). 

6.2 Low-mass leptoquark and squark limits 

To illustrate the sensitivity of this search, 95% CL. upper limits on the cross section 
times the branching ratio, <J/3e q , for F = and F = 2 leptoquarks are shown in Fig. 6; for 
the e + p case, only the subsample (65 pb" 1 ) with the higher y/s of 318 GeV is used. Upper 
limits on \ eqi *J [3 m are shown in Figs. 7 and 8 for F = and \F\ — 2 scalar and vector 
LQs, assuming resonantly produced leptoquarks as described by the BRW model. Since, 
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for sufficiently large LQ masses, the cross section is dominated by electron valence-quark 
fusion, only e + p (e~p) data were used to determine F — (\F\ — 2) LQ production limits. 
Similar considerations hold for the results shown for the e — r channel in Figs. 9 and 10. 
For couplings with electromagnetic strength (A e?1 = \i qfj = 0.3 ~ y/4ira), LQs with 
masses up to 299 GeV are excluded (see Tables 1 and 2). Alternatively, for a fixed Mlq 
of 250 GeV, values of \ eqi (3 m and of \ eqi (3 Tq down to 0.010 and 0.013, respectively, 
are excluded (see Tables 3 and 4). 

Constraints on \ eqi ^J~f3i~ q for Sy 2 and for Sq can be interpreted as limits on A' ljl y / /3^Z^ 
and ^'iiks/ Pd k ->eq f° r & ano - d k -R-Parity- violating squarks of generation j and k, respec- 
tively [39]. 

6.3 High-mass leptoquark and squark limits 

Tables 5 and 6 show the 95% C.L. limits on X eqa A^/Mlq (third row of each cell) for 
F = and \F\ = 2 high-mass leptoquarks coupling to eq a and \iqp. Limits were evaluated 
for all combinations of quark generations a, (3, except when a coupling to a t quark is 
involved. Tables 7 and 8 show the corresponding limits for LQs coupling to eq a and rq^. 

Limits for Sy 2 LQs can also be interpreted as limits on Ay a A^/M~ for a w-type squark of 
generation j, where % = 2, 3 is the generation of the final-state lepton (/x or r). Similarly, 
limits for Sq LQs can also be interpreted as limits on A' lafc A^ fc /Mj for a <i-type squark of 
generation k. 



7 Comparison with limits from other experiments 
7.1 Low-energy experiments 

There are many constraints from low-energy experiments on lepton-flavor-violating pro- 
cesses coming from muon scattering and rare lepton or mesons decays [8]. Most can be 
converted into limits on \ eqa A^/Mlq for massive scalar or vector leptoquark exchange. 
In Tables 5-8, the limits from such measurements are compared to the constraints from 
this analysis. For the e — \i transition, such indirect limits are very stringent and ZEUS 
limits are better only in a few cases involving the c-quark. In the e — r channel, ZEUS 
improves on the existing limits for many initial- and final-state quark combinations, espe- 
cially when a quark of the second or third generation is involved. Assuming X eqi = Xe qp , 
ZEUS limits on low-mass LQs can be compared to the limits from low-energy experiments. 
In Figs. 7 and 8, limits on X eqi as a function of the LQ mass are compared to the limits 
from e — n conversion in nuclei and from rare K- and B- meson decays. ZEUS limits are 
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better or competitive with indirect limits up to ~ 250 GeV when the quark in the final 
state is of the third generation. In Figs. 9 and 10, the corresponding limits for the r case 
are shown compared to constraints from rare r, B or K decays. ZEUS limits improve on 
low-energy results in most cases. 



7.2 LFV and leptoquark searches at colliders 

Tevatron limits are complementary to those from HERA since the cross sections at pp 
colliders do not depend on the Yukawa coupling, and LQs are assumed to couple only 
with one lepton generation. Therefore, such experiments are sensitive to only a subset 
of the interactions considered here. The CDF and D0 collaborations exclude scalar LQs 
coupling exclusively to fiq with masses up to 202 GeV [40] and 200 GeV [41], respectively. 
CDF performed an analysis searching for leptoquarks which couple exclusively to the 
third generation of leptons and excluded LQs with Mlq < 99 GeV if /3 r & = 1. The D0 
collaboration looking for vvbb final states excluded LQs with masses below 94 GeV if 
j3 u b = 1. The CDF collaboration also performed a search for a narrow resonance decaying 
to two charged leptons of different generation [42], observing no deviation from the SM 
expectation. 

Searches for LFV interactions, not mediated by LQs, were performed by LEP experiments, 
looking for e/i, er and fir production in e + e~ annihilation at the Z° peak [43]; the OPAL 
collaboration extended the search to higher energy using LEP2 data [44]. Also in this 
case, no significant deviation from the SM expectation was found. 



8 Conclusions 

The data taken by the ZEUS experiment at HERA in e + p and e~p interactions at center- 
of-mass energies of 300 GeV and 318 GeV during the years 1994-2000 corresponding to an 
integrated luminosity of 130 pb~ x were analyzed for lepton-flavor violation. Searches in 
both fj, and r channels were performed. No evidence of lepton-flavor- violating interactions 
was found. For masses lower than the center-of-mass energy, limits at 95% C.L. were set 
on \ e q 1 \[~fiiq for leptoquark bosons as a function of the mass. For a coupling constant 
of electromagnetic strength (X eqi = \i qf} = 0.3), mass limits between 257 and 299 GeV 
were set, depending on the LQ type. For Mlq = 250 GeV, upper limits on \ eqi */ (3 m 
{Xeqx-^Wq) in the range 0.010-0.12 (0.013-0.15) were set. 

For LQs with M L q 3> y/s, upper limits on X eqa \ mj3 / and X eqcc X Tq g/MlQ were calcu- 
lated for all combinations of initial- and final-state quark generations. 
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Some of the limits also apply to R- Parity- violating squarks. In many cases, especially in 
the r-channel, ZEUS limits are more stringent than any other limit published to date. 
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LQ type 


CL 
°l/2 


CL 
°l/2 


qR 
°l/2 


V L 


yR 


V R 




/i-channel limit on Mlq( GeV) 


273 


293 


293 


274 


278 


296 


299 


r-channel limit on Mlq ( GeV) 


270 


291 


291 


271 


276 


294 


298 



Table 1: 95% CL. lower limits on Mlq for F = LQs in the fi- and the 
r-channels assuming X eqi = Xg q0 = 0.3. 



LQ type 


oL 
D 


qR 
°0 


cR 

D o 


st 


V L 

V l/2 


yR 
V l/2 


yR 


/i-channel limit on Mlq ( GeV) 


278 


284 


261 


281 


269 


289 


289 


T-channel limit on Mlq( GeV) 


275 


281 


257 


278 


265 


287 


286 



Table 2: 95% CL. lower limits on M L q for \F\ = 2 LQs in the \i- and the 
r-channels assuming X eqi = Xi qf} = 0.3. 



LQ type 


CL 
°l/2 


CL 
°l/2 


qR 
°l/2 


v L /v* 


yR 


Vi 


/i-channel limit on X eqi \J fi m 


0.054 


0.021 


0.019 


0.037 


0.015 


0.010 


r-channel limit on X eqi \V {3 Tg 


0.066 


0.026 


0.024 


0.046 


0.019 


0.013 



Table 3: 95% CL. upper limits on A egi y7^g f or F = LQs with mass Mlq = 
250 GeV in the fi- and the r-channels. 



LQ type 


oL I qR 


qR 
°o 


st 


V L 

V l/2 


yR 
V l/2 


% 


/i-channel limit on X eqi \Jf3 w 


0.047 


0.12 


0.041 


0.080 


0.030 


0.033 


r-channel limit on X eqi (3 Tq 


0.058 


0.15 


0.049 


0.10 


0.038 


0.042 



Table 4: 95% CL. upper limits on X eqx ^f^ q for \F\ = 2 LQs with mass Mlq = 
250 GeV in the fi- and the r-channels. 
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e ->• n ZEUS e±p 94-00 F = 


a/3 


°l/2 

e~u 
e+u 


cR 
°l/2 

e (u + a) 

e+(u + d) 


°l/2 

e~d 
e+d 


V L 

e~d 
e+d 


yR 

e~ d 
e+d 


yR 

e~u 
e+u 


e (v 2m + a) 
e+{V2u + d) 


1 i 


fj,N -» eTV 
5.2 x 1CT 5 
1.2 


uN -> eAT 
2.6 x 1CT 5 
1.0 


/xTV eAT 
5.2 x 10~ 5 
1.7 


fiN -» 
2.6 x 10~ 5 
1.0 


fiN -> eAf 
2.6 x 10~ 5 
1.0 


fiN e/V" 
2.6 x 10~ 5 
0.8 


fiN -> eAT 
0.8 x 10~ 5 
0.4 


1 2 


D —y fxe 
2.4 


K —y fie 
2 x 1(T 5 
1.0 


AT — ► /ie 
2 x 10~ 5 
1.8 


A' — > /ie 
1 x 10~ 5 
1.2 


AT — > /ie 
1 x 10~ 5 
1.2 


D — > fie 
1.2 


A' — > /xe 
1 x 10~ 5 
0.5 


1.3 


1.0 


1 3 


* 


B ^ /ie 
0.4 
1.8 


B -> /ie 
0.4 
1.9 


B — > fie 
0.2 
1.5 


B ^ fie 
0.2 
1.5 


* 


B — > /ie 
0.2 
1.5 


2 1 


D — > /ie 
2.4 
3.6 


K ^ fie 
2 x 10~ 5 
2.4 


K — > /te 
2 x 10~ 5 
3.1 


K — ► /ie 
1 x 10~ 5 
1.3 


if — > /te 
1 x 10~ 5 
1.3 


D ^ fie 
1.2 
1.2 


K —> fie 
1 x 10~ 5 
0.6 


2 2 


//AT -► eTV 
9.2 x 10~ 4 
5.7 


//AT -> eiV 
1.3 x 10~ 3 
3.1 


jLtTV -> eN 
3 x 10~ 3 
3.8 


/iAf -► eAT 
1.5 x 10~ 3 
1.9 


fiN ^ eN 
1.5 x 10~ 3 
1.9 


/zA" -> e/V" 
4.6 x 10~ 4 
2.8 


/xA" -> eAT 
2.7 x 10~ 4 
1.1 


2 3 


* 


S -> peif 
0.3 
4.3 


B -> ueAT 
0.3 
4.2 


B — > /2eAT 
0.15 
2.9 


S -> /xeA' 
0.15 
2.9 


* 


£ -> /ieAT 
0.15 
2.9 


3 1 


* 


B ^ fie 
0.4 
4.4 


B —> fie 
0.4 
4.4 


Kb 
0.12 
1.5 


B ^ fie 
0.2 
1.5 


* 


V ub 
0.12 
1.5 


3 2 


* 


B -> ueAT 
0.3 
5.8 


-B — > /leA' 
0.3 
5.8 


S -> fleK 
0.15 
2.2 


S -> ueAT 
0.15 
2.2 


* 


B —> fieK 
0.15 
2.2 


3 3 




fiN ^ eN 
1.3 x 10~ 3 
7.6 


/xTV -> eAT 
3 x 10" 3 
7.6 


fiN -> eN 
1.5 x 10~ 3 
3.9 


fiN ^ eN 
1.5 x 10~ 3 
3.9 


* 


fiN eAT 
2.7 x 10~ 4 
3.9 



Table 5: Lzmzfo at 95% C.L. on e *"/ 9/3 for F = LQs, m units of TeV~ 2 . 

LQ 

The first column indicates the quark generations coupling to LQ — e and LQ — \x, 
respectively. ZEUS results are reported in the third line (bold) of each cell. The 
low- energy process providing the most stringent constraint and the corresponding 
limit are shown in the first and second lines. The ZEUS limits are enclosed in a 
box if they are better than the low-energy constraints. The cases marked with * 
correspond to processes where the coupling to a t quark is involved. 
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Table 6: Lzrmfe at 95% C.L. on eg °' 190 for F = 2 LQs, in units of TeV~ 2 . 

LQ 

The first column indicates the quark generations coupling to LQ — e and LQ — p, 
respectively. ZEUS results are reported in the third line (bold) of each cell. The 
low- energy process providing the most stringent constraint and the corresponding 
limit are shown in the first and second lines. The ZEUS limits are enclosed in a 
box if they are better than the low-energy constraints. The cases marked with * 
correspond to processes where the coupling to at quark is involved. 
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Table 7: Limits at 95% C.L. on eq " ^ for F = LQs, m units of TeV~ 2 . 

1V1 LQ 

The first column indicates the quark generations coupling to LQ — e and LQ — r, 
respectively. ZEUS results are reported in the third line (bold) of each cell. The 
low- energy process providing the most stringent constraint and the corresponding 
limit are shown in the first and second lines. The ZEUS limits are enclosed in a 
box if they are better than the low-energy constraints. The cases marked with * 
correspond to processes where the coupling to a t quark is involved. 
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Table 8: Limits at 95% C.L. on ZV"' 3 for F = 2 LQs, in units of TeV~ 2 . 

LQ 

The first column indicates the quark generations coupling to LQ — e and LQ — t, 
respectively. ZEUS results are reported in the third line (bold) of each cell. The 
low- energy process providing the most stringent constraint and the corresponding 
limit are shown in the first and second lines. The ZEUS limits are enclosed in a 
box if they are better than the low-energy constraints. The cases marked with * 
correspond to processes where the coupling to at quark is involved. 
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Figure 1: (a) s-channel and (b) u-channel diagrams contributing to LFV processes. 
The subscripts a and (3 denote the quark generations, and I is either a \i or a r . 



21 



ZEUS 



w 

-4—* 

c 

> 

11110 



(a) 



b i « gjj i , , , , i , , , , i , , , , i , 




• ZEUS e±p 94-00 

Background MC 
. LFV n M LQ =240GeV 



10 20 30 40 50 60 70 80 90 100 

E-P z + A^ (GeV) 



20 40 60 80 100 120 

f x (GeV) 



to 

4-» 

CD 
> 

11110 



(c) 

't 



2 3 4 5 6 7 8 9 10 11 12 

1/2 




^ t /\fi;(GeV) 



20 40 60 80 100 120 140 160 



Figure 2: Comparison between data (dots) and SM MC (solid line): (a) E — Pz + 
A M; (b) ft, (c) ft/y/El and (d) polar angle of the muon, 9^, after the ^-channel 
preselection. The dashed line represents the LFV signal due to a scalar LQ, with 
Mlq = 240 GeV , with an arbitrary normalization. 



22 




Figure 3: Comparison between data (dots) and SM MC (solid line) for the 
variables used in the r discriminant: (a) — log(-R mcan ); (b) — \og(R Tma ); (c) 
-log(l -L mean ); (d) -log(L rms ); (e) number of subjets, N suhi ; (f) jet mass, M jet; 
after the r -channel preselection (hadronic r decays). The dashed line represents 
the LFV signal with arbitrary normalization. 
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Figure 4: Distribution of (a) the discriminant, T>, and (b) A<p, after hadronic r 
decay 'preselection. The dots represent the data while the solid line is the SM predic- 
tion from MC. The LFV signal distribution for two different LQ masses, 240 GeV 
(dashed line) and 600 GeV (dash-dotted line), are also shown with arbitrary nor- 
malization. The distribution of A<j) for the Mlq = 600 GeV LQ, which is similar 
to the Mlq = 240 GeV LQ A0 distribution, is omitted. The leptonic decay of the 
tau, or the tau jet outside the CTD acceptance, leads to events with A0 > 160°. 
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Figure 5: A0 distribution of the events with T> > 0.9 after hadronic r decay 
preselection. Dots represent data while the solid line is the SM prediction from 
MC. The dashed line represents the signal with arbitrary normalization. The small 
fraction of the signal (~ 5%) with A<p > 160° is due to events that have the jet from 
the t outside the CTD acceptance. The two events from data that have A0 = 72° 
and A<p = 126° are the two events found in a previous ZEUS search for isolated t 
lepton events [31]. 
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Figure 6: The 95% C.L. upper limits for a f3 \ q as a function o/Mlq for scalar (full 
line) and vector (dashed line) LQs: (a) F = 2 LQ — > \iq; (b) F — LQ — > fiq; (c) 
F = 2 LQ — > rq; (d) F = LQ — > rq. A subset of e + p data (99-00, corresponding 
to the higher center- of -mass energy, 318 GeV) has been used to obtain figures (b) 
and (d). 
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Figure 7: Limits for F = low-mass LQs in the [i channel obtained from e + p 
collisions. The upper plots show 95% C.L. limits on \ eqi x W~(3^q for (a) scalar and 
(b) vector LQs. In the lower plots, ZEUS limits on X eqi for a representative (c) 
scalar and (d) vector LQ are compared to the indirect constraints from low-energy 
experiments [8], assuming X eqi = X m/3 . 
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Figure 8: Limits for F = 2 low-mass LQs in the /i channel obtained from e~p 
collisions. The upper plots show 95% C.L. limits on \ eqi x W~j3^q for (a) scalar and 
(b) vector LQs. In the lower plots, ZEUS limits on X eqi for a representative (c) 
scalar and (d) vector LQ are compared to the indirect constraints from low-energy 
experiments [8], assuming X eqi = \ Ml3 . 
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Figure 9: Limits for F = low-mass LQs in the r channel obtained from e + p 
collisions. The upper plots show 95% C.L. limits on \ eqi x -J [3 Tq for (a) scalar and 
(b) vector LQs. In the lower plots, ZEUS limits on \ eqi for a representative (c) 
scalar and (d) vector LQ are compared to the indirect constraints from low-energy 
experiments [8], assuming X eqi = \ T q p - 
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Figure 10: Limits for F = 2 low-mass LQs in the r channel obtained from e~p 
collisions. The upper plots show 95% C.L. limits on \ eqi x j3 rq for (a) scalar and 
(b) vector LQs. In the lower plots, ZEUS limits on X eqi for a representative (c) 
scalar and (d) vector LQ are compared to the indirect constraints from low-energy 
experiments [8], assuming X eqi = X Tqp . 
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